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The discovery of microwave propagation well beyond the radio horizon 
in the late forties created an immediate interest in its communication possi- 
bilities. The first decade (1950 - I960) was characterized by a vigorous 
experimental and theoretical program aimed at mapping out those propagation 
features which were most relevant to the design of radio relay systems. Long 
term averages of mean signal level and their rough scaling with frequency and 
distance were of primary interest. Highly reliable signal levels were empha- 
sized at the expense of unusually large, though infrequent signals. The fine 
structure of the signal was explored in a rough way, in order to set conserva- 
tive bounds on the communication channel capacity. Space correlations of the 
signal were measured so as to establish appropriate separation distances for 
diversity transceivers. Theoretical explanations were tied either to a layered 
reflection model or to turbulence scattering schemes. The single scattering 
far field or cross section approximation was used almost exclusively to describe 
the latter, while signal variability had to be embroidered onto the stable layer 
theories by qualitative concepts of moving glint or reflection points. The two 
theories were first considered as competitive explanations of the same phe- 
nomenon, although gradually it was appreciated that a variable mixture of the 
two was probably responsible. Whatever the deficiencies of this program, one 
must admit that the combination of rough experiment and theory met the 
communication engineering needs of the first decade rather well. 

The second decade beginning in I960 is probably best described as a 
renaissance of the scientific interest which first stimulated the subject. The 
attractive possibility of using short term radio measurements as a means for 
inferring the instantaneous structure of the atmosphere through which the 
waves had passed was gradually appreciated. Important new experimental 
techniques also became available to improve such capabilities. Very stable 
oscillators were developed which made it possible to measure round trip 
signal phase on scatter paths (1). The time domain correlation or RAKE 
technique used by the Sylvania group(^) provides a unique means for examining 
particular range slices of the scattering volume. Multiple arrays of receiving 
antennas were used by the Stanford group ( 3 ) to study the azimuth and elevation 
fine structure of the propagation volume by the electronic synthesis of very 
narrow beams. Over and above these measurement techniques, data processing 
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capabilities have improved importantly since the early days when much of the 
data was hand reduced. The capacity to do high speed digital processing of 
wide band signals is central to the availability of substantial amounts of pre- 
cision data on signal fine structure. Taken together, these experimental 
techniques provide a means for studying the atmospheric structure in three 
space dimensions and time with exciting precision. Nevertheless, our ability 
to design specific measurement 'programs to exploit this capability depends in 
large measure on our intrinsic understanding of the phenomenon. 

Unfortunately, our theoretical understanding of the propagation mechanism 
has not kept pace with experimental progress. No adequate description of the 
signal fine structure has been developed thus far. We must still depend on 
two qualitative, competitive theories (layers vs. scatter) to explain the general 
features of the transmission, as discussed in Dr. Cox's review ( 4 ). Either a 
reconciliation or fusion of the two theories is needed to suggest the next step 
in an evolving understanding of the propagation. It may even be that the two 
models are simply different ways of describing the same basic mechanism 
which is as yet undiscovered. It is vital to separate the geometrical consider- 
ations of the propagation path from the physics of the layers and/or turbulent 
irregularities. Gaussian correlation models, chosen for their integrability , 
are no longer an adequate starting point. Nor is it appropriate to try to force 
all available data into universal turbulence theories which may or may not be 
primarily responsible at various times. Finally, the familiar restrictions on 
stationarity , homogeneity, and isotropy should be relaxed in a consistent way. 

There has been a great deal of advanced theoretical research on line-of- 
sight propagation, both in Russia and the United States, stimulated in large 
measure by optical propagation using lasers. However, the geometries and 
signal statistical considerations for line -of-sight and transhorizon propagation 
are sufficiently different to make the recent research largely irrelevant to 
scatter propagation. 

The general features of a theory for single scattering by turbulent 
irregularities which meets most of the above needs were described in 1959 
The fundamental physical notion is that the received signal is the integrated 
result of scattering/reflection events throughout the common volume. This 
leads naturally to an integral (equation) relation between the measured signal 
and the refractive structure producing it. 

ts(e') - yv GCv.-r'i A? Cf, A (l) 

V 

where GfR^) is the free space Green's function, V the common volume, A ^ 
the electromagnetic wavelength, R the receiver and T the transmitter locations. 
Formally, this is just the Born approximation for scattering by a dielectric 
perturbation AiJ. . The problem is that is the function that we wish to 
infer, not assume. Furthermore, it is usually a stochastic function of position 
and time, which is defined only by its statistical averages. The essential 
trick is to introduce a three dimensional Fourier transformation 
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in which the turbulence wavenumber ]<. decomposition of the turbulent 
irregularities Ck.4^ is related to the more familiar wavenumber spectrum 
^ by: 

< W,v>W,^> -= 3E3 £ u\ - (3) 

CZTsV 5 

If we substitute (2) into (1) and interchange the orders of integration, we 
achieve the crucial separation of propagation geometry and turbulence physics, 

^ \ \ ^ ^ V <2- CCe v O ( 4 ) 

V 

Physics Geometry 


Finally one can establish an explicit expression for the average scattered 
power by taking the magnitude of (4) and introducing Equation (3). 

\) cf c e EoC'T.-r') I (5) 
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It is this expression, rather than the cross section approximation, which 
provides the starting points for a more precise theory. 

All of the electromagnetic and path geometrical features of the propaga- 
tion are contained in the bracketed quantity in (5) , which appears as a weight- 
ing function or kernel of the turbulent spectrum S ( K )• The frequency and 
distance dependence of the scattered field are implicit in this weighting function 
via their appearance in the incident wave E 0 , the Green's function, and their 
convolution with <2jty>vK’'C r over the scattering volume V. One can also exploit 
Equation (5) to describe more complicated experiments. For instance, by 
inserting appropriate Delta functions in the volume integration, it is possible 
to isolate those elements of the received signal corresponding to particular 
multipath components. An inhomogenous field in which the intensity of 
turbulent irregularities varies with altitude is relevant to the real atmosphere 
and can be included in the above expression as height dependent shaping 
factors ( te enc? ) in the geometrical integrals. The influence of a 

transmitter antenna pattern is contained naturally in the incident wave E 0 , 
and can be included for the receiver by multiplying the free space Green's 
function G(R, r) by the pattern function describing ray propagation from the 
various scattering points r to the receiver R. The usual approximation for 
narrow beam geometries, in which the integrations are replaced by V times 
an average value of the integrand, can thus be checked. By further 
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complicating the geometry of Equation (5), one can describe beam swinging 
within or off the great circle plane. Unfortunately, this formalism has not yet 
been reduced to explicit expressions useful for comparison with experiment, 
and these calculations represent an ambitious program. However the important 
point to make is that it is essentially an exercise in integral calculus, one 
which can reliably provide the relationship between the scattered power and 
the still arbitrary spectrum S(.<) for any propagation path. The hope, of course, 
is that such integral relations can be inverted to deduce the spectrum from 
the measured quantities. 

It is of some interest to note that the bracketed weighting function in 
Equation (5) is a function of the three vector components of the wavenumber 
vector This occurs because of the different ways in which the propagation 

geometry integrals emphasizes the x, y, and z projections of K. via the 
scattering mechanism. Through its argument, this produces a directional 
emphasis of the spectrum S( k„ \t^), which must be considered together with 
non-isotropy of the turbulent fiela itself. The possible confluence of flat 
layers with large horizontal anisotropic blobs having short vertical correlation 
has been suggested before ( 6 ) ( and these two observations may provide an 
avenue to a common understanding of the phenomenon. 

The basic signal statistical distribution of the received field should also 
provide a mechanism for distinguishing between: (1) scatter, and (2) layers 

plus scatter or glint. The usual assumption that the amplitude is Rayleigh 
distributed is not confirmed in any detail by experiment - nor should it be 
expected. The central limit theorem for the contributions of a large number 
of scattering elements does tell one that the orthogonal signal components of 
the scattered field 


G-s - XCt ^ + (6) 

ought to be distributed as follows: 

exo -\ * V \ (7) 

^ 2 v? <r z Y i. - ^ 

where is the cross correlation coefficient of x and y and CT is the 
variance of each. Converting to polar amplitude and phase coordinates and 
integrating over the phase 
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If the correlation jo were zero, as usually assumed, this would reduce to 
the Rayleigh distribution. However, one can generalize Equation (^1 to 
write an explicit expression for 
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In a mature theory this should be calculated as a function of geometry and 
compared with experiment. However, we must also note that a constant 
vector A added to a scattered component x + iy also produces a non- Rayleigh 
distribution. Starting from (6) and setting /° = o for contrast and temporary 
simplicity, we find 

G'cA - <10 Tj\ no) 

which is what one ought to expect from a steady layer reflection plus a 
turbulent scatter component or a random glint from the layer. The hope is 
that a detailed comparison of short sample experimental data with the two 
predictions can distinguish between the mechanisms. 

The relationship of signal amplitude and phase to one another at succeed- 
ing instants of time provides one of the most useful measures of signal fine 
structure. Included in such relationships are the time auto -correlation of 
amplitude and phase, plus their corresponding power spectra. The derivatives 
of these quantities lead to signal fading rates and doppler shifts, both of which 
are now available experimentally. The joint probability distribution of two 
time-displaced complex signals 

EG (. t') be. , -v ' *0^ ^ 

and 

x t ^ (ll) 

has been given before in its most general form (9). By converting to polar 
signal coordinates, it is possible to write an explicit expression for the joint 
distribution of amplitude Ri, R 2 , and phase and <fi^. From this one can 
compute all relevant experimental quantities by integration; viz <C , 

<((, , etc. The coefficients in this probability distributions are 
functions of propagation geometry and the time displacement C . Explicit 
expressions for the coefficients which are time-displaced generalizations of 
equations like Equation (8) have been established ( 5 ) using the appropriate 
generalization of Eq. (3). 
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Here Lc is the drift velocity of a frozen irregular structure and 
describes the time correlation of the turbulent internal rearrangement of the 
structure. Using this general decomposition, it is possible to write 




x CGc.'cN 
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and so on for the other components of the moment matrix which defines the 
coefficients in the general probability distribution function (5). Notice that 
these expressions also provide explicit separation of the turbulent physics and 
the path geometry. In point of fact, the integrals which occur in the time- 
displaced signal component correlations are the same as those required for 
the ordinary variances of x and y plus the cross correlation • 

The novel physical ingredient is the time autocorrelation of the turbulent 
structure C (K ,•£) describing the self motion. Most qualitative theories of 
fading on scatter paths have ignored this and dealt only with the horizontal 
wind vector LL . The intuitive justification is that the vertical transport 
velocity is negligible and that only off great circle scattering elements have 
significant horizontal components of ~\£ via the weighting integrals in (13). 

This may be, but it is worth verifying by explicit calculation. The self-motion 
has no such projection, so that its vertical component is effective in producing 
time variability in the great circle plane. There is some reason (5) to believe 
that C (K,-£, ) ought to be a function of , which would suggest a fading 

rate variation with carrier frequency of , as compared to a linear 

variation for straight drift motion fL . Experiment shows a sliding variability 
between the two extremes, suggesting a mixture. This kind of experiment 
ought to provide an exceedingly good means for studying the effect and relative 
strengths of the various propagation mechanisms. Furthermore, time vari- 
ability measurements are probably easier to perform than beam swinging 
experiments because they do not require long times for antenna alignment and 
give one a chance to examine true snapshots of the atmospheric structure. 

A final word about space and frequency correlations is in order. The 
basic probability distribution for signal components gathered at displaced 
antennas or at the same site on separated frequencies is given by the same 
expression referenced (5) above for time-correlated signals. The coefficients 
which occur in such distributions change by virtue of the modification of the 
geometrical integrals viz: 
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while the separated frequency expressions involve different electromagnetic 
wavenumbers in the Green's function and incident field. The calculation of 
these coefficients is part of the same unified analytical program suggested 
above. 

Having pointed the way to a wider theory, it is appropriate to emphasize its 
deficiencies. Firstly, it refers only to the turbulent scatter component and 
not to a steady layer reflection. A comparable sophisticated theory is needed 
for the latter, plus a consistent means for predicting its time variability, 
space correlation, etc. We have also assumed single scattering in the 
common volume, and some have suggested that multiple scattering may be the 
dominant mode . While our expressions do allow for a non-isotropic 
turbulent spectrum Vt) depending individually on the three components of 
l<r, we have no clear means for dealing with a non- stationary process. 
However, the primary deficiency of this proposal is that it is a plan and not a 
finished thesis. Perhaps the availability of new experimental data and the 
renaissance of genuine scientific interest in the subject will stimulate such 
research. 
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